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ABSTRACT 



Relativistic binary pulsars, such as B1534+12 and B1913+16 are character- 
ized by having close orbits with a binary separation of ~ 3 Kq. The progenitor 
of such a system is a neutron star, helium star binary. The helium star, with a 
strong stellar wind, is able to spin up its compact companion via accretion. The 
, ^ ! neutron star's magnetic field is then lowered to observed values of about ~ 10^° 

Gauss. As the pulsar lifetime is inversely proportional to its magnetic field, the 
possibility of observing such a system is, thus, enhanced by this type of evolution. 
O ! We will show that a nascent (Crab-like) pulsar in such a system can, through 

accretion-braking torques (i.e. the "propeller effect") and wind-induced spin-up 
"r^ i rates, reach equilibrium periods that are close to observed values. Such processes 

occur within the relatively short helium star lifetimes. Additionally, we find 
O I that the final outcome of such evolutionary scenarios depends strongly on initial 

^ [ parameters, particularly the initial binary separation and helium star mass. It 

^ ' is, indeed, determined that the majority of such systems end up in the pulsar 

"graveyard", and only a small fraction are strongly recycled. This fact might 
^ I help to reconcile theoretically expected birth rates with limited observations of 

5^ \ relativistic binary pulsars. 
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1. INTRODUCTION 

Two of the four known High Mass Binary Pulsar systems, (HMBP's) - PSR B1913+16 
(Hulse & Taylor 1975) and PSR B1534+12 (Wolszczan 1990) - have short orbital periods 
(~ 10 hours). Such systems, upon their eventual mergers, are considered to be important 
sources of gravitational wave radiation which may be measured by the next generation of 
detectors. As a result, it becomes desirable to understand the evolutionary processes that 
lead to their formation. Note that we do not consider PSR B2127+11C (Prince et al. 1991) 
which resides in the globular cluster, M15: globular cluster sources may have completely 
different evolution mechanisms which are dominated by dynamical interactions. 

The link between relativistic binaries and their original 0/B main- sequence progenitor 
systems are beheved to be wide High Mass X-Ray Binaries (or Be/HMXB's). The standard 

evolutionary scenario following the X-Ray phase (Bhattacharya & Van den Heuvel 1991; 
Van den Heuvel & Van Paradijs 1993) predicts that the neutron star enters the hydrogen 
envelope of its giant companion. Common Envelope (CE) evolution ensues as the compact 
object spirals in, creating dynamical friction and ultimately expelling the envelope. The 
orbit is then tightened, leaving a helium star, neutron star binary. However, Chevalier 
(1993) showed that a neutron star in CE evolution would likely form a black hole. Brown 
(1995a) confirmed this scenario, showing that hypercritical accretion forces > IMq onto the 
neutron star, sufficient to form a black hole. 

Brown's alternate scenario for the formation of relativistic binaries involves, instead, a 
double helium star binary (Brown 1995a; Wettig & Brown 1996). If the progenitor 0/B 
supergiants (ZAMS mass ~ 24Mq) are initially very close in mass (within four percent), the 
two stars will burn helium at the same time. Thus it is possible for the neutron star to avoid 
moving through the envelope of the secondary. Although CE evolution takes place, it does 
so with two helium stars. Furthermore, a natural explanation is given as to why the pulsar, 
which gains mass by accretion, is the heavier star in the binary. Such a result is supported 
by observations. Notice that either scenario leads to a neutron star, helium star binary as an 
intermediate stage. Thus, the evolutionary scenario discussed here is generally applicable. 

The pulsar born into such a binary is Crab-like with a strong magnetic field 3 — 10 x 10^^ 
Gauss, and a short spin period, 30-50 ms. Small orbital separations (1-3 R0) and strong 
helium star winds ensure heavy accretion onto the neutron star, causing its magnetic field 
to lower two orders of magnitude and spinning it up further. Competing with wind-fed 
accretion is the so-called "propeller effect." This mechanism exerts a spin-down torque on 
rapidly rotating neutron stars as material is prevented from accreting. We shall see that 
such a scenario is, indeed, consistent with observations of relativistic binaries such as PSR 
1913-1-16. Note that diminished magnetic field strengths lengthen the observable hfetime of 
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a pulsar (r = P^/2S^). Thus, an "observability premium" is given to recycled pulsars. 

In §2, we discuss our model of the evolution of a hehum star and neutron star in a close 
binary. Particularly, we examine how the spin period and magnetic field of the neutron star 
as well as the orbital separation of the system change with time. Since many parameters 
are involved, an analytic solution is not readily available. Instead, in §3, we discuss results 
of a computer code which was set up to analyze the model. Additionally, since initial 
conditions are not well understood, we also discuss the effect of variation of parameters on 
possible outcomes. It is determined that the final outcome of the binary evolution strongly 
depends on such parameters (particularly initial orbital separation and helium star mass) 
and, furthermore, the number of systems which would lead to an observable, recycled pulsar 
is highly constrained. 



It has been shown (Langer 1989; Woosley et al. 1993) that when a helium star loses 
matter from an enhanced stellar wind, the mass-loss rate is dependent on its total mass. 
Furthermore, it is assumed that when very massive helium stars explode in a supernova, there 
is a tendency to dissociate the binary. As we are solely interested in modeling progenitors 
to relativistic binaries, we consider, as a first approximation, the lower-bound mass limits 
from Woosley et al. (1993), i.e. M^e = -5 x lO-^Mgyr-^ (MHe/M0)^•^ 

However, it has recently been argued (Brown et al. 2001) that this rate is too high by 
a factor of 2-3. This is supported by polarization measurements of Thomson scattering in 
helium stars as well as the observed scaling of mass-loss rates with orbital period changes 
(St.-Louis et al. 1993; Moffat & Robert 1994). Using polarization measurements of V444 
Cygni, (Mne = 9.3Mq) St.-Louis et al. (1993) found Mhc = 0.75 x lO-^Meyr"!. (Note that 
they use a slightly different value for the terminal wind velocity. See eq. [5]) Extrapolating 
this rate to all helium star masses and assuming the same scaling as in Woosley et al., we 
therefore find: 



The helium star nuclear burning lifetime is also dependent on its total mass. In an 
approximation to evolutionary calculations made by Paczyhski (1971) and Habets (1986), 
Pols et al. (1991) estimate a functional form for helium star lifetimes as 
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DESCRIPTION OF THE MODEL 
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1.6 < M/Mq < 4.8 
M/Mq > 4.8 



(2) 



Habets (1986) has shown that hehum stars with M < 2.2Mq become white dwarfs and, 
thus, are not considered here. Also note that, as we assume the progenitor to our model to 
be a double helium star binary, we only use half of the total helium star lifetime, given by 
equation (2), in our calculations. 

We then look to the question of how much of the ejected helium star matter is actually 
accreted onto its compact companion (Mx). The captured mass rate, /c, is defined as the 
fraction of mass captured by the neutron star's gravitational field and is given by /c — 
—Mca,p/M-iie where Mcap is the mass capture rate. Note that, as Mne is negative, fc is always 
greater than or equal to zero. Similarly, the accreted mass fraction, /a is given by the fraction 
of mass transferred that is actually accreted onto the surface of the compact star. Here, /a 
= Mx/Mcap. As matter can not be accreted faster than the Eddington limit, we define /a = 
1 if Mcap < MEdd but /a = -/^Edd/ -^cap for super-Eddiugtou transfers. In sum, we see: 



Alternatively, one may define the parameter a as the total fraction of matter that is lost 
from the system. In that case, a — 1 — fcfa.- 

Assuming a standard Keplerian orbit with binary separation, a, the neutron star moves 
relative to the helium star with a circular orbital velocity 



In an approximation to numerical work done by Habets (1986) we estimate that the helium 
star wind velocity at the position of the neutron star obeys the relation given by. 



where the helium star radius is given by Rub/Rq — 0.22(MHe/Mo)°-^ and Voo ~ 2000 km 
s~^. Generally, the helium stellar wind will move out radially, orthogonal to v^, so that the 
pulsar experiences the wind moving at a relative velocity of Vr — {v^ + v'^Y^'^. 

Adopting the standard accretion mechanism (Bondi 1952), we assume incoming matter 
gets captured near the so-called accretion radius, Rg. This is defined as the point where the 



Mx = -/c/aMHe 



(3) 




(4) 




(5) 
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wind velocity (relative to the neutron star) is equal to the escape velocity of infalling matter. 
Thus, 

i?g = — ^ (6) 

i/j. 

It should be noted that Bondi (1952) defines the accretion radius as Rg = 2GMx/ {v^ + cf), 
where Cg represents the speed of sound in the plasma. However, for the parameters of this 
model, this may be considered negligible and ignored. 

Finally, we can use geometry to estimate the captured mass fraction, /c. We determine 
it to be the fraction of a sphere of radius, a, which occupies the area enclosed by the accretion 
radius of the neutron star, i.e. /c = TrR'^/iira^. With equation (6), this becomes: 

/c = ^'"^^ , (7) 
a2 (^2 + vlf 



2.1. Orbital Evolution 

In order to determine how the orbital separation changes with time, one needs to model 
how angular momentum is transferred in the binary. Neglecting spin angular momentum 
and assuming circular Keplerian orbits one finds Jorb = A*a^<^ where u) = ^ G{M^ + Mne) / 
and /X is the reduced mass. Recalling that Mx — {a — l)MB.e and differentiating, one can 
easily show 

a _ 2 Jorb _ 2MHe 
a Jorb 

Since matter is being lost from the system, it is clear that total angular momentum 
will not be conserved. It is then necessary to make some assumptions about how angular 
momentum is transferred from one star to another. Reasonably, one could expect that the 
actual mass being passed from the helium star to the neutron star (i.e. fc dM-ne) is transferred 
conservatively. We then assume that the fraction (1 — fc)dMiie will leave the system with the 
specific angular momentum of the helium star, jHe- Thus dJ; = (1 — fc) (iMneine = (1 ~ fc) 
dM]ie{MyL/ ^toifo?^- Similarly, we expect that /c(l — fa)dMiic leaves the system with the 
specific angular momentum of the neutron star, dJf = /c(l — /a) (iMHc(MHc/Mtot)^a^u;. 
Together, these assumptions yield a differential equation for the orbital angular momentum 
of the system, Jorb = Ji + Jf , 



Mne a f Mne \ 
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Hc 



■/orb ^ M] 
Jorb + MHe 



MHe 



l-/c)T7^+/c(l-/a)- 



Ml 



He 



(9) 



Defining the mass ratio, q = M^/Mffe, and combining equations (8) and (9), we find 



a _ q[q-fc{2-2q' + f,q)] 
a [q + fcf.W + Q) 



(10) 



2.2. Emitter Phase 



Initially, the newborn pulsar is characterized by a rapid spin (P= 30 — 50 ms.) and a 
strong dipole magnetic field {Bg ~ 3 — 5 x 10^^ G). Such enhanced dipole radiation pressure 
might be sufficient to keep the wind plasma from being accreted onto the neutron star. In 
this phase the compact star will behave like an isolated radio pulsar, spinning down with 
time as rotational kinetic energy is converted to dipole radiation energy. Such systems have 
been extensively studied (Gunn & Ostriker 1969; Goldreich & Julian 1969) and we model 
dipole radiation pressure for a neutron star of radius and spin period, fl as 



'''^^''^ ~ ~ ^^^^ 

The stopping radius, Rs, is defined as the point where dipole radiation pressure is suf- 
ficient to balance wind pressure from the helium star (Urpin et al. 1997). If the stopping 
radius is less than the accretion radius, then the neutron star will behave like an isolated 
emitter, spinning down via the Gunn-Ostriker mechanism. However, if R^ > Rg, then it is 
possible for the pulsar to accrete. We estimate wind pressure to be given by ~ PwV^ 
where pw, the plasma density, is approximately given by pw ~ |-^He|/47ra^Vw Finally, setting 
— -Prad(-Rs)) wc get an expression for the stopping radius 



R^ = A i 12 

^'6c4|MHe|^;w 



Note that as no accretion may take place at this point (/c = 
field will remain constant (see eq. [18]). The spin period, P = 
to the standard relation (Gunn & Ostriker 1969) 



/a = 0), the pulsar magnetic 
27r/Q will increase according 
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Here, I = k^ME? represents the moment of inertia of the neutron star with k"^ ~ 0.4 (Prakash 
et ah 2000). 



2.3. The Propeller Phase 

According to equation (12), the stopping radius is proportional to the inverse square of 
the period. Thus, over a relatively short amount of time, the pulsar will spin down sufficiently 
such that accreting matter may interact with the magnetosphere of the neutron star. Recall 
this occurs at the point where the stopping radius falls inside the accretion radius. 

We assume, for simplicity, that the magnetosphcric boundary of the neutron star is 
defined where the ram pressure of infalling matter is balanced by the neutron star's magnetic 
dipole pressure (Lamb, Pethick & Pines 1973) Assuming spherical inflow (but see Ghosh & 
Lamb, 1979a,b), the magnetospheric radius is, thus: 



= I ^ (14) 

Once matter couples to the neutron star's magnetic field, the interaction's effect on 
the overall spin evolution depends on the balance between centrifugal and gravitational 
accelerations. If the pulsar rotation is, initially, too fast, the neutron star will eject infalling 
plasma, propelling it away tangentially while simultaneously losing angular momentum in 
the process. This "propeller" effect (Pringle & Rees 1972; Illarionov & Sunyaev 1975) can 
be parameterized by a fastness parameter, = i7/f2k(-Rm) where Q\^{R^) = y^GM^R^a, 
is the Kcplcrian angular velocity at the magnetospheric boundary. Clearly, for a fastness 
parameter greater than unity, the propeller mechanism is initiated. 

The resulting angular momentum loss (and associated increase in spin period) can be 
represented by a propeller torque, N^^.^^ = IQ = -Mcap-Rmt'esc(-Rm) = -Mcapy/2GM^R^. 
Consequently, the spin-rate decreases according to the relation: 

^ = ^V^gK^ (15) 
It is instructive to note that it is also possible to use energy methods to determine the 
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propeller torque. Over time, rotational kinetic energy of the neutron star will be transmitted 
through shocks to the wind plasma falling near the magnetospheric boundary (Fabian 1975). 
Consequently, this gas will heat up and be dispersed when it attains escape velocity {V ~ 
Ksc = ^j2dMjR^). Thus we find, E = md = -l/2M,apKsc = -GM^M,,JR^. Hence, 
7VpE„p ^m= -GM^M,,jR^n = -[GM^M,,jRA{Rm)]uJ-\ 

Upon examination of the two possible propeller torques, we see 



<op = --^VgMA. = (16) 

Thus, for extremely fast rotators, it is much more difficult for an energy propeller to change 
the spin period of the neutron star significantly. The reality of the magnetospheric interaction 
is very complex, so it is hard to say which mechanism is inherently more accurate. Therefore, 
we use both the energy and angular momentum propellers and compare results. 



2.4. Accretion 

For fastness parameter less than unity, co-rotating matter is able to accrete to the surface 
of the neutron star. At this point, as the helium star wind carries angular momentum, we 
expect the pulsar spin period to decrease over time. Thus, the pulsar is said to be recycled. 

It is not initially clear whether spherical or disk accretion will dominate over the course of 
the helium star burning time. Shapiro & Teukolsky (1983) argue that if the intrinsic angular 
momentum per unit mass of accreted gas, jaco exceeds the specific angular momentum of an 
element in a circular Keplerian orbit near the magnetospheric radius, iKcp(-Rm) = y/GM^Rm, 
then disk accretion will dominate. Otherwise, we may treat the mass-transfer as being 
(nearly) spherical. Thus, the necessary prerequisite for a disk to form is jacc ^ ixep where it 
may be shown (Shapiro & Teukolsky 1983), jacc — (l/2a)i'x-Rg and V:^. and Rg are given by 
equations (4) and (6) respectively. 

It then becomes desirable to include the effects of magnetic torques within the disk on 
the overall spin-rate (viscous torques may be ignored here). This complex problem was first 
discussed by Ghosh & Lamb (1979a,b) where they found that for slow rotators {ujs <S 1) 
magnetic coupling may enhance spin-up torques by as much as forty percent. For cus i, 

the opposite is true and magnetic effects might actually oppose the spin-up. Following 
Ghosh & Lamb, we define, for disk accretion, A^^cc — n{uJs)Myj[R^ where the dimensionless 
coefficient, valid for < < 0.9, is given by 
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. . 1.39 [4.03 (l-^,)"-'"'- 0.878]} , , 

n{u,) ^ !^ — ^ 17 

Next, we discuss the evolution of the neutron star's magnetic field. Recent observations 
and analyses seem to strongly indicate that mass accretion in binary systems is directly 
correlated with magnetic field decay in neutron stars (but see Wijers, 1997) . Although there 
exist many possible accounts of physical mechanisms that would explain this phenomenon 
(Konar & Bhattacharya 1997), here we rely on the empirical model of Shibazaki et al. (1989) 
and leave a more detailed analysis of accretion driven magnetic field decay for later work. For 
typical initial values of dipole magnetic field strengths, i.e. B12 ~ 1 — 10, it was determined 
that one can make the following empirical approximation to magnetic field evolution: 



m = — (18) 

Here, Ms ~ 12.5 x 10^*^ Mq is a typical scaling parameter. This parameter has been set 
to qualitatively agree with the observed magnetic fields of Low Mass X-Ray Binaries. Over 
their lifetimes, LMXB's can accrete up to ~ O.IM0 and have typical magnetic fields on the 
order of 5 x 10^ Gauss. Assuming an initial field strength of ~ 5 x 10^^ G, we see that 
equation (18) roughly gives the desired result. 



2.5. Valving and Equilibrium 

As we have shown, the line of demarcation between spin-up (propeller effect) and spin- 
down (accretion) phases of evolution is parameterized by the fastness parameter, lJs- For 
a very fast rotator, a;^ ^ 1 and the spin period increases with time as the magnetic field 
remains constant. However, as the neutron star spins down, matter co-rotating with the 
neutron star at the so-called co-rotation radius, Rc, will spin down as well and Rc oc P^/^. 
As the propeller phase continues, the co-rotation radius will increase until it is finally greater 
than the magnetopsheric radius, given by equation (14). At this point, the neutron star will 
allow co-rotating matter to fall to its surface and the accretion phase begins. 

However, once accreting, the magnetic field begins to decay. By equations (14) and (18), 
we see that R^ is, roughly, a monotonically decreasing function of time and as accretion 
continues, the magnetospheric radius diminishes, possibly allowing the propeller phase to 
resume. 

In sum, we expect a kind of oscillation or "valving" to occur between accretion and 
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propeller phases of the neutron star's evolution, until equilibrium is restored. For spherical 
mass transfer at equilibrium the co-rotation radius should coincide with the magnetospheric 
boundary. In this case, one can estimate the equilibrium spin period of the neutron star: 



So, for a canonical neutron star accreting spherically at the Eddington limit, Pgq = 
(13.6 ms)B^Q^. As the magnetic field is brought down with each cycle of accretion, Peq is 
brought down as well. 

The equilibrium point for disk accretion is generally more complicated. Upon examina- 
tion of the Ghosh and Lamb function, (eq. [17]), we find that n{ujs) is undefined for uJs = 1. 
In this case, we define the critical ratio, ujc = 0.5050. Since n{u!c) ~ —1 and for all > ujc 
magnetic torques are sufficient to force the neutron star to spin down, we define cUg = uj^ as 
the turnover from a propeller to an accretor. I.e., we define all n{uis > cue) — —1. 

In conclusion, the overall effect of the evolution, over the total helium burning time, is 
expected to result in (a) lowering the neutron star magnetic field, (b) spinning the pulsar 
up to or near the millisecond range, and (c) a slight widening of the orbit. 



As discussed in the introduction, and expanded on in the previous section, there are 
several parameters that need to be accounted for in order to evolve the binary system with 
time. We are, thus, led to several coupled differential equations for which there is no analytic 
solution. Instead, we now discuss the results of a computer code set up to analyze the model. 

At each time step within the total helium burning time (eq. [2]), the stopping radius, 
accretion radius, magnetospheric radius and fastness parameter were calculated in order to 
determine the predominant phase of the evolution (i.e. emitter, propeller or accretor). The 
differential equations of §2 were then solved numerically, using a simple Euler scheme, and 
the calculated parameters were appropriately updated. 

Initial conditions such as stellar masses, neutron star spin period, magnetic field strength 
and orbital separation for relativistic binary progenitors are largely unknown so the simula- 
tion was run for a wide range of values. Specifically, the initial helium star mass was varied 
in the range 2.2Mq < Mne < 15M0. Helium stars with masses below the lower hmit tend 
to evolve into white dwarfs (Habets 1986) whereas stars with masses far above 15M0, will 




(19) 
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RESULTS OF THE NUMERICAL CALCULATIONS 
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dissociate the binary upon supernova. Neither possibihty would lead to relativistic binary 
systems. It has been shown (Brown et al. 2001) that due to their strong stellar wind, naked 
helium stars in this mass range generally end their lives as neutron stars and not as black 
holes. Similarly, as we are interested in the progenitors of close binaries, the initial orbital 
separation was chosen to have the (somewhat arbitrary) maximum value of IORq. 

The criterion determining the minimum orbital separation was defined such that the 
helium star stay within its Roche lobe, i.e. i?L < Rrc- If a helium star were allowed to 
expand beyond its Roche lobe, unstable mass-transfer would commence. Such a scenario 
would generally not lead to the formation of a HMBP and is, thus, excluded. Egglcton 
(1983) has shown that for a mass-ratio, q, (eq. [10]) the Roche-lobe radius of a star is given 

by 



i?L _ 0-19 ^ ^ ...^ 

a ~ 0.6 + ?2/3in(l + 5-1/3) -^^^^ 



By equation (20), we see that flmin = Ruc/fiq) and is, therefore, a function of helium 
star mass. Since Rne/Re = 0.22{Mne/MQ f-^ , a,^,^{2.2MQ) = 0.8435/?© and a^U'^^MQ) = 
1.9URq. 



The initial properties of the nascent neutron star were constrained to a somewhat smaller 
range with the Crab pulsar (PSR 0531-21) taken as a prototype for our model. Canonical 
values were given for neutron star masses, radii and moments of inertia (1.471/©, 10^ cm,~ 10^^ 
g-cm^ respectively). Initial spin periods ranged from 30 - 50 ms. and initial dipole field 
strengths varied on the order 1 — 10 x 10^^ Gauss. 

Several plots of the time evolution of spin period and magnetic field under various 
conditions follow and certain general trends are evident. A sharp rise in spin period with 
time (with no corresponding change in magnetic field) indicates that the primary factor 
dominating the evolution is the propeller effect. However, when the spin period lengthens 
to some critical value, the propeller mechanism gives way to accretion. At this point, the 
period and magnetic field strength both decrease with time. 

For sufficiently large initial binary separations and/or low helium star masses, accretion 
may not occur at all during the total helium burning time. In such a case, the dipole field of 
the neutron star remains constant and its overall spin period actually increases with time, 
due to the propeller effect and electromagnetic spindown. Many times, this spin-down is 
sufficiently large to put the neutron star in the pulsar graveyard. 

Figure 1 examines the time evolution of the pulsar spin period and magnetic field for 
both the energy and angular momentum propellers (see §2.3) under typical conditions. Inital 
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parameters are given by Mue,! — 4.OM0, Bi — 5.0 x lO^^Gauss, Pi — 50.0 ms, and ai = 
1.50it!o- The total helium burning time is calculated, using equation (2), to be 6.246 x 10^ 
years. 

Clearly, the final results and overall evolution are sensitively dependent on the type of 
propeller mechanism used. For the case of the angular momentum propeller, the neutron 
star quickly enters the accretion phase at t — 2.42 x 10^ years. Consequently, it ends up 
with a short final period (99.74 ms) and low magnetic field strength (4.169 x 10^° G). Thus, 
the given initial conditions, coupled with the J-propeller mechanism, is sufficient to recycle 
the pulsar. 

Things are quite different for the energy propeller, however. Prom Figure 1, we see the 
neutron star never accretes at all and the magnetic field remains constant throughout the 
evolution. Weak torquing inhibits the rapid spin-down seen for the J-propeller. Instead, the 
period remains roughly constant with Pf = 121.1 ms. Additionally, we see that, although 
the angular momentum propeller reaches a period close to the equilibrium period (eq. [19]) 
early in its evolution (~ 4 x 10^ years), the E-propeller never does. 

For both cases, the helium star loses 0.4848 Mq by a steady wind and, as anticipated, 
most of this mass is lost. The neutron star accretes 1.494 x 10^'^ for the J-propeller 
while mass transfer is completely non-conservative for the E-propeller. Note that for this 
low-mass helium star, mass-capture is always sub-Eddington and /a = 1 for all time. Finally, 
we see that the type of propeller mechanism has httle effect on the widening of the orbit 
(of = 1.644it!o, 1.645i?0 respective). 

Figure 2 illustrates the case where a low-mass hehum star (SM©) is in a relatively wide 
orbit (oi = 8Rq). The initial properties of the neutron star are otherwise identical to the 
previous case. Here we see that accretion does not set in until i ~ 5.7 x 10^ years for the 
angular momentum propeller while the energy propeller, again, does not accrete at all over 
the total evolution time of 9.897 x 10^ years. 

For the first phase of the evolution, the weak wind of the hehum star is insufficient to 
push past the stopping radius of the neutron star. The pulsar acts like an emitter, spinning 
down electromagnetically. Then at t ~ 4.9 x 10^ years, the concavity of the curve in the P 
vs. T diagram changes. At this point the neutron star has spun down sufficiently for the 
helium star wind to begin to interact with its magnetosphere. However, the pulsar is still a 
fast rotator and the propeller phase commences. 

The type of propeller mechanism determines the next stage of the binary's evolution. 
For the case of the angular momentum propeller, the neutron star rapidly spins down to a 
maximum pulse period greater than 43 seconds. At t = 5.65 x 10^ years, accretion finally 
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begins but equilibrium is never achieved. Only 9.880 x lO^^M© is actually accreted to the 
surface (AMne = — O.3636M0) and the final spin period and magnetic field strength of the 
neutron star are 33.14 seconds and 2.798 x 10^^ Gauss respective. Clearly the pulsar has 
entered the graveyard and this system would not be observed as a relativistic binary. 

Once again, for the energy propeller, weak torquing prevents significant changes in 
the spin period and, although the neutron star never accretes at all (magnetic field remains 
unchanged), the final spin period is 298.9 ms. As a result, the neutron star can be observed as 
a radio pulsar (yet it is clearly not recycled). For both scenarios, the final orbital separation 
is 8.72Oi?0. 

Finally, we examine the case where a massive helium star lies in a close orbit with a 
neutron star. In this case we expect there to be heavy accretion and recycling. In particular, 
consider a 12Mq helium star in orbit with a canonical neutron star (Pj = 50 ms, 5i = 5 x 10^^ 
G) with an initial binary separation of 3.Oi?0. 

Figure 3 illustrates the resulting evolution of spin and magnetic fields for the angular 
momentum propeller. We sec there is, indeed, heavy accretion and after a sharp propeller 
cycle, equilibrium is achieved at P{ = 60.41 ms. Additionally, the neutron star accretes 
2.086 X 1O~^M0 which is sufficient to bring the field down to 2.993 x 10^° Gauss in a time 
of 2.668 x 10^ years. The final orbital separation is 3.833i?0. 

Mass transfer is spherical throughout the evolution and, due to the high initial mass 
of the helium star, wind loss is supcr-Eddington until t = 2.54 x 10^ years. The minimum 
period of 59.35 ms occurs near t ~ 2.1 x 10^ years where steady accretion gives way to 
valving. At this point, the period slowly increases and closely tracks the (time-dependent) 
equilibrium period. 

The energy propeller again prohibits accretion onto the neutron star. Through electro- 
magnetic spindown and propeller torque, the pulsar gradually increases its spin to a final 
period of 83.71 ms. All of the hehum star wind matter (2.944M0) is lost while the magnetic 
field remains constant. 

For sufficiently heavy accretion, it is possible that slight changes in initial orbital sepa- 
ration may result in significant differences in the neutron star's spin period evolution. Con- 
sider the case of a lOM© helium star in orbit about a neutron star with initial magnetic field 
7.5 X 10^^ G, and initial spin period 35 ms. Assume an angular momentum propeller effect 
only. Figure 4 shows the resulting time evolution for oi = 1.957?© and 2.05Rq respectively. 
We find that the binary in the closer initial orbit actually has a longer final spin period. 
Such a result seems counterintuitive and needs to be examined more closely. 
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When Qi — 2.O5i?0, accretion is spherical throughout the mass transfer phase, lasting 
from t = 6.1 X 10^ years until the time of the helium star's supernova ed, t — 2.98 x 10^ years. 
The final, equilibrium spin period is 48.4 ms. However, for the case when aj = 1.95i?0, 
accretion changes from disk-type to spherical. The short, initial, electromagnetic spin-down 
phase is nearly identical for both cases, as is the initial propeller phase. Then, at t ~ 5.6 x 10^ 
yr, spherical accretion begins. This continues, uninterrupted until t ~ 4.2 x 10"^ years. At 
this point, a second, smaller, propeller phase begins and continues until t — 4.25 x 10^ 
years. Valving and disk accretion sets in at this time, allowing the spin period to gradually 
decrease with time. Finally, at t — 2.46 x 10^ years, jacc falls below jKep and spherical 
accretion commences again. This results in a much faster drop in spin-period with time. 
The neutron star spherically accretes until the helium star supernova and the final spin 
period for this case is 55.0 ms. 



4. DISCUSSION 

Next we examine the results of a contour plot which parameterizes the final spin period 
and magnetic field as functions of initial orbital separation and helium star mass. Standard 
initial values were again chosen for the neutron star with Mx,i = 1.4M0, Bi = 5x 10^^ G, and 
Pi = 50ms. As discussed in §3, the initial helium star mass range is 2.2M0 < M-n^i < 15M0 
whereas the initial orbital separation varied in the range flmin < fli < 1O7?0. Here, amin is 
defined to be R^e/ f{q) (eq. [20]). The output of the contour plots representing the angular 
momentum propeller and energy propeller appear in figures five and six, respectively. 

From the slope of the contours in Figure 5a, it is clear that, for close orbits, the final 
spin period most strongly depends on the initial orbital separation, and less on the helium 
star mass. However, for wider (initial) orbits, both parameters play an important role on 
the final outcome. 

If neutron stars act as angular momentum propellers, then for any orbit with initial 
separation greater than ~ 6-R0, the pulsar will not only be unrccyclcd but it will sit in the 
graveyard as well. Therefore, there is a strong constraint placed on progenitors of relativistic 
binary pulsars such as 1913-1-16. 

Recall that pulsar spin-down times are proportional to P^/B"^. Thus, one expects an 
"observability premium" for low-field, fast rotators. I.e., due to longer lifetimes outside 
the graveyard, one should expect to observe a higher proportion of such systems. Indeed 
this is true of the two relativistic binaries (three including 2127-l-llC). It seems that, upon 
examination of figure 5a, that in order to bring the field below 5 x 10^° Gauss (O.OlSi), 
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the initial binary separation must be a maximum of ~ ARq. Additionally, if one wishes to 
observe a system with spin periods on the order of 50 ms, even closer orbits are necessary, 
with ao ~ 2 — 3Rq. And this is only true for the most massive of helium stars, i.e. those with 
A^Hc,i > 8Mq. Thus for angular momentum type propellers, we conclude that only a minute 
portion of the overall parameter space represented in figure four will produce relativistic 
binary pulsars such as PSR's 1913+16 and 1534+12. 

The parameter space of final outcomes is much more highly constrained for the energy 
propeller (figures 6a,6b) than for the angular momentum propeller. It seems, from examining 
the spin period contour plot in figure 6a that it is extremely difficult (if not impossible) for 
a neutron star to spin down into the graveyard by the E-propeller mechanism. Even in the 
weak- wind, wide orbit limit, Pf never rises much higher than ~ 0.5 seconds. At the other 
extreme, we find that it is also very difficult to force the neutron star to accrete at all for 
this scenario. Figure 6b shows that mass transfer, parameterized by magnetic field decay, 
occupies only a tiny portion of M-a space. 

Prom equation (16), we see that the propeller torque is inversely proportional to the 
fastness parameter. Thus, unlike the J-propeller, fast rotators are unable to easily change 
their spin period. Thus, one may conclude that the energy propeller is highly inefficient at 
recycling pulsars. 

This characteristic becomes a significant factor in determining whether an E-propeller 
will accrete at all or not. As with the J-propeller, slightly changing the initial conditions 
(particularly the initial orbital separation) may make a profound difference in the final out- 
come of the pulsar's evolution (although, for different reasons). The result of this condition 
is the dense band of contours near the Omin line in figure 6a. 

For a specific example of why this occurs, consider a 6Mq helium star orbiting a neutron 
star with initial parameters Bi = 4x 10^^ G and Pi = 70 ms. Figure 7a shows a plot of final 
spin period as a function of initial separation in the range amin = 1.27P0 < Oi < 2.5Rq. A 
sharp peak exists at the point, Oi = 1.567^©. Here, Pf reaches a maximum value of 3.063 
seconds. For > 1.6Pq, Pf falls dramatically to a level of about ~ 85 ms at a; = 2.1Pq. 
At this point there is a much more gradual increase in final spin period with increasing 
separation. 

What causes the sudden shift in final spin period? As figure 7b shows, a slight change in 
initial orbital separation will have an effect on the fastness parameter and, thus, the overall 
evolution of the propeller phase. Here, we examine the time evolution of the spin period for 
the three different initial orbital separations - 1.525^©, 1.565^©, and 1.6O5P0. As the initial 
orbital separation increases, the propeller torque will decrease by a small amount and, as a 



-16- 



result, it takes longer for the neutron star to reach its maximum period. For ai — 1.525Rq, 
the spin period reaches a maximum of 3.09 seconds at t = 3.43 x 10^ years. Afterwards, the 
accretion phase begins and continues until the end of the helium star lifetime at 4.044 x 10^ 
y. where Pf = 177.7 ms. A relatively long accretion phase allows the magnetic field to decay 
to a final value of 1.358 x 10^^ G. So, we see a slightly recycled neutron star for the given 
initial conditions. 

Changing the initial separation to 1.565i?0 increases the fastness parameter and, con- 
sequently, decreases the propeller torque. Now the peak docs not occur until nearly the end 
of the evolution (t = 4.03 x 10^ years). As there is little time for the accretion phase, the 
magnetic field does not change much (i?f = 2.6 x 10^^ G) and the final spin period is 2.33 
seconds. Thus, the pulsar winds up in the graveyard and is not observable. Finally, by 
increasing the initial separation to 1.6O5i?0 we can completely eliminate the accretion phase 
altogether. Here, the pulsar peak never occurs and the neutron star can not be recycled. 
The magnetic field remains constant throughout the evolution (Pf = 234.9 ms). 

Next we consider the case of PSR J1518+4904 (Nice et al, 1996). Of all four known 
HMBP's, this one is the most recycled (P = 40.94 ms, logP = 9.1). However, it resides in a 
relatively wide binary with an orbital period of 8.634 days (a ~ 25^©). Unfortunately, our 
model can not account for such properties. 

From figure 5 a (J-propeller), we see that if ai > 6Rq, the pulsar is destined to end 
up in the graveyard. Even for the most massive helium star companions, the neutron star 
will never accrete and, instead, a large propeller torque will sharply increase its spin period. 
Such a trend diminishes for ^ IORq since, for very wide binaries. Pg > Rg for most (if 
not all) of the evolution. The helium star wind is too weak to interact with the neutron 
star's magnetosphere at such separations, and the pulsar spins down by the Gunn-Ostriker 
mechanism. Thus, we would expect that J1518+4904 would not be observable nor would it 
undergo magnetic field decay. 

Although the evolutionary history of PSR 1518+4904 is unknown, it is clear that our 
model does not account for all possibilities. We, therefore, must assume its evolution was 
different than the other HMBP's. One possible evolutionary scenario which may account for 
the observed properties of 1518+4904 is reverse case C mass-transfer (Kippenhahn & Weigert 
1967). As low-mass helium stars tend to have significantly extended envelopes (Habets 1986), 
even an initially wide binary may undergo sufficient mass-transfer to initiate recycling during 
the helium shell-burning stage. Assuming the neutron star survived the resulting spiral-in 
phase, it may be possible to end up with a recycled pulsar in a wide binary. Thus, wc find 
that two possible improvements to our current model would be an accurate treatment of 
case C mass-transfer and a more detailed physical model of magnetic field decay. 
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5. CONCLUSION 

Relativistic binary pulsars, such as B 1534+ 12 and B 19 13+ 16 arc characterized by hav- 
ing close orbits (a --^ 3Rq), with recycled pulse periods and magnetic fields (P ~ 30 — 60 ms, 
log 5 ~ 10). A single case of a wide HMBP (B1518+4904) exists. We do not consider PSR 
2127+llC as it resides in the globular cluster, M15 and may have a different evolutionary 
history. 

We assume wind- fed mass transfer is responsible for recycling the observable neutron 
star in a HMBP. In lieu of a more detailed physical mechanism, we use the empirical model 
of Shibazaki et al. (1989) to account for magnetic field decay. Such a model assumes changes 
in the pulsar magnetic field are proportional to the amount of of mass accreted to its surface. 

We find that for close initial orbits (ai ~ 2 — 3 Rq) and an initial helium star mass in 
the range SMq < MHe,i ^ 15Mq, we are able to reproduce the observed spin periods, orbital 
separations, and magnetic fields for HMBP's 1913+16 and 1534+12. Brown et al. (2001) 
find that helium stars in this mass range will generally end their lives as neutron stars. 
Whereas such high mass He stars are relatively rare, they do not undergo a significant red 
giant phase during helium shell burning. Such a phase would, in fact, make it more difficult 
to form a binary pulsar for the less massive systems because reverse case C mass-transfer 
would generally lead to the formation of a black hole (Fryer & Kalogera 1997). Finally, we 
note that with such a high initial helium star mass, a kick velocity upon supernova would be 
required if the system is to remain bound. This is consistent with earlier calculations made 
by Burrows & Woosley (1986). 

Our model was unable to account for the properties of PSR 1518+4904 and we speculate 
that another mechanism such as reverse case C mass-transfer is responsible for its properties. 

Cycles of accretion coupled with the propeller effect allow the neutron star to come 
to be recycled in a time consistent with helium star nuclear lifetimes. Of the two possible 
propeller mechanisms proposed, only the angular momentum propeller is efficient at recycling 
the HMBP progenitors to observed properties. 

As pulsar lifetimes are inversely proportional to their magnetic fields (r ~ P'^/B^), re- 
cycling and magnetic field-decay lengthen the observable lifetime for a neutron star. Thus, 
an "observability premium" is introduced for relativistic binary pulsars. It is, furthermore, 
discovered that the final outcome of the binary evolution strongly depends on initial con- 
ditions (particularly MHe,i and Oi). Thus, the number of possible observable systems are 
constrained to a small region of the overall parameter space of initial conditions. 
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Fig. 1. — Time evolution of the pulsar spin period (solid lines) and magnetic field (dashed lines) for both 
the angular momentum and energy propellers. For both cases, Si = 5 x 10^^ G, Pi = 50 ms, and Oi = l.SOi?©. 
Initial helium star mass is 4.0Mq. For the .J-propcUcr. wc sec evidence of recycling with Pf = 99.74 ms and 
Bi = 4.169 X 10^" Gauss. The neutron star accretes 1.494 x IQ~^Mq over an evolution time of 6.246 x 10^ 
years. For the E-propeller, weak torquing inhibits a rapid spin-down and, consequently, there is no accretion. 
The final spin period is 121.1 ms while the magnetic field does not change. 




Fig. 2. — Time evolution for the case of a low-mass helium star (SMq) in a wide orbit (oi = SP©). Other 

initial conditions arc the same as in figure 1. For the J-propeller, the pulsar will end up in the graveyard 
with a final spin period of 33.1 seconds. Only 9.88 x IQ'^Mq is accreted onto the surface of the neutron 
star (AMne = -0.3636Mq) and there is little field decay (Pf = 2.798 x 10^=^ Gauss). For the E-propeller, 
once again, there is no accretion as evidenced by the constant magnetic field. Here, the final spin period is 
298.9 ms and the neutron star is observable as a radio pulsar. For both scenarios, the total evolution time 
is 9.897 X 10^ years and the orbit widens until Of = 8.72OP0. 
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Fig. 3. — Spin and magnetic field evolution for the case ai = 3i?o,MHe,i = 12Mq (angular momentum 
propeller only). All other initial conditions are the same as in previous cases. There is heavy accretion 
and, after a sharp propeller cycle, the neutron star spins up to an equilibrium period of 60.41 ms. The field 
strength is, subsequently, lowered to a final value of 2.993 x 10^° Gauss in a time of 2.668 x 10^ years. The 
final orbital separation is S.SSSi?©. For the E-propeller (not shown), there is no recycling and Pf = 83.71 
ms. 
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Fig. 4. — Here we see the time evolution of the spin period for a canonical neutron star (Pi = 35 ms, 

Bi = 7.5 X 10"'^^ G ) in a close orbit with a lOAf© helium star (J- propeller). For Oj ~ 2.05Rq. there is 
continuous spherical accretion from t — 6.1 x 10^ years until t[ = 2.98 x 10^ years. The final equilibrium 
spin period is 48.4 ms. However, by slightly changing the initial conditions, such that = l.QSP©, the time 
evolution becomes quite different. Here, spherical accretion dominates from 5.6 x 10^ yr. until 4.2 x 10^ yr. 
This leads to a second (smaller) propeller phase and a cycle of disk accretion. Finally, at t = 2.46 x 10^ 
years, accretion again becomes spherical and Pf = 55 ms. 
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Fig. 5. — Contour plots of (a) final spin period and (b) final magnetic field strengths as a function of initial 

separation and helium star mass for the angular momentum propeller effect. Initial values are given using 
the Crab pulsar (PSR 0531-21) as a prototype. Therefore, we have set Bi = 5 x 10^^ Gauss and Pi = 50 ms. 
Variation of these parameters do not change the final outcome much. The initial helium star mass varied 
in the range 2.2 < Mue < 15Mq and initial orbital separation range was amin < a < IORq. The amin line 
is defined by auiin = Ruc/Hq) (sec text). The lower right portion of the graph indicates heavy recycling 
whereas the upper left contours represent neutron stars that have spun into the graveyard. 
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Fig. 6. — Contour plots of (a) final spin period and (b) final magnetic field as functions of initial orbital 

separation and helium star mass for the c;nergy propeller. Initial conditions and parameter space are the 
same as for the angular momentum propeller. Note that final values are much more strongly constrained 
and, it seems, it is very difficult for the E-propeller mechanism to either spin the neutron star into the 
graveyard or sufficiently recycle it to form a relativistic binary such as PSR 1913+16. In figure (6a) we see 
that final spin periods near the flmin line are sensitively dependent on the initial conditions {see text for an 
explanation). From (6b) we find that accretion, and, subsequent magnetic field-decay is absent for the vast 
majority of the parameter space in the case of the E-propeller. 
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Fig. 7. — (a) A plot of final spin period as a function of initial separation in the range 1.27 < Oi < 2.5Rq 
(Energy propeller, Bi = 4 x 10^^ G, Pi = 70 ms, Muc.i = QMq). A sharp maximum occurs at the point 
tti = 1. 567 Rq where Pmax = 3.06 seconds. For ai > 1.6i?0, Pf falls dramatically to about ^ 85 ms at 
tti = 2.1i?0. At this point, the final period gradually rises again with increasing separation, (b) A plot of 
the time evolution of the spin period for three different initial orbital separations - 1.525i?0, 1.565i?0, and 
1.6O5i?0. See text for details. 



